Magnetic Feshbach resonances in collisions of closed-shell X S molecules 
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Magnetic Feshbach resonances play a central role in experimental research of atomic gases at 
ultracold temperatures. A major thrust of current research is to create an ultracold gas of diatomic 
alkali-metal molecules in the ground rovibrational state of the ground electronic X E state. Can 
ultracold 1 T, molecules be controlled by means of magnetic Feshbach resonances? Unlike alkali 
metal atoms, 1 T, diatomic molecules have no unpaired electrons. The response of X E molecules to an 
external magnetic field is determined by the spin structure of the atomic nuclei. We present quantum 
scattering calculations illustrating that nuclear spin relaxation and Feshbach resonances in ultracold 
molecule-molecule collisions for 1 E molecules in a magnetic field are highly probable. We show that 
certain Zeeman states of X E molecules exhibit broad resonances and weak inelastic scattering, which 
makes these states ideal for experiments exploiting magnetic field control of ultracold molecular 
gases. 

PACS numbers: 



Much of the success of experimental research with ul- 
tracold atomic gases is due to the possibility of tuning 
the scattering length of ultracold alkali metal atoms by 
means of magnetic Feshbach resonances [![. The scat- 
tering length near a Feshbach resonance undergoes rapid 
variation as a function of the magnetic field, which has 
been used for experimental studies of Bose-Einstein con- 
densation (BEC) bosonic superfluidity quan- 
tum magnetism |5[, many-body spin dynamics [6[, Efi- 
mov states 0, Bardeen-Cooper-Schrieffer (BCS) super- 
fluidity H and the BEC-BCS crossover 0, E3|. A ma- 
jor thrust of current experimental work is to extend 
this work to ultracold molecules 
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Of particular in- 
terest arc polar alkali metal dimcrs in the ground 1 Y, 
electronic state produced by photoassociation of ultra- 
cold atoms. Several experiments have recently demon- 
strated the creation of ultracold alkali metal dimers KRb 
12| in the ro- vibrational ground state of the electronic 
gro und state and rapid progress is being made with RbCs 

iEE|. 

Polar molecules offer long-range dipolar interactions, 
which can be used for new exciting applications such as 
the study of dipolar crystals [IH and ultracold controlled 
chemistry [li, 17 1. While the long-range interactions be- 
tween polar molecules can be controlled by an external 
electric field [3, [HI], many applications, such as quantum 
simulation of spin-lattice models [20[ or the study of dipo- 
lar effects on BEC and the BEC - BCS cross-over [211 ] . 
require independent control over short-range and long- 
range interactions. The short-range interactions could, 
in principle, be controlled by magnetic Feshbach reso- 
nances. However, unlike alkali metal atoms, 1 E molecules 
have no unpaired electrons and it is currently believed 
that the binary interactions between 1 Y, molecules can- 
not be controlled by a magnetic field. 

For two interaction partners A and B, a magnetic Fes- 
hbach resonance occurs when a bound state of the AB 



complex has a different magnetic moment from that of 
the isolated particles so that the energy of the bound 
state can be tuned across an A + B collision threshold by 
varying an external magnetic field. The width of the res- 
onance is proportional to the coupling between the bound 
and continuum states. For collisions of two alkali-metal 
atoms in the ground electronic state 2 S, these couplings 
arc determined by the energy difference between the sin- 
glet ( 1 E) and triplet ( 3 E) electronic states of the collision 
complex and strong hyperfine interactions between the 
unpaired electron spin and the nuclear spin [l|. Because 
1 S molecules have a closed electronic shell, the interac- 
tion properties of two 1 H molecules are determined by 
a single potential energy surface. The response of ^ 
molecules to an external magnetic field is determined by 
the magnetic moments of the nuclei. The nuclear spin 
states of atoms or molecules with closed electronic shells 
are generally well protected from the effects of collision 
dynamics. However, we show in the present work that 
a subtle interplay of couplings induced by the molecule 
- molecule interaction potential and the interactions in- 
duced by molecular rotations gives rise to magnetic Fes- 
hbach resonances that can be used to control ultracold 
collisions of 1 E molecules. 

Figure 1 shows the Zeeman energy levels of the 
87 Rb 133 Cs(X 1 E + ) molecule in the ground rotational 
state. At zero magnetic field, the energy levels are la- 
beled by the quantum number F of the total angular 
momentum including the nuclear spins. At high mag- 
netic fields, the energy levels can be labeled by the quan- 
tum numbers M^b and Mq s specifying the projections 
of the nuclear spins of Rb and Cs on the magnetic field 
axis. Recent experiments demonstrated that ultracold 
alkali metal dimers can be selectively prepared in any 
Zeeman level of the ground rotational state 22j. In 



the present work, we examine the strength of collision- 
induced couplings between the Zeeman states shown in 
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Figure 1. These couplings lead to nuclear spin relax- 
ation and magnetic Feshbach resonances in molecule - 
molecule collisions. As follows from the equations pre- 
sented below, if two 1 £ molecules were initially prepared 
and remained during a collision in the rotational ground 
state with N = 0, the nuclear spin states would have 
been completely decoupled from translational motion of 
the molecules. However, the molecule - molecule interac- 
tion potential induces couplings between rotational states 
that modify the hyperfine interactions during collisions. 
These couplings may give rise to magnetic Feshbach reso- 
nances and induce the nuclear spin relaxation in molecule 
- molecule collisions. 

In order to calculate the energy levels shown in Figure 
1, we diagonalize the following Hamiltonian 
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Magnetic Field (Gauss) 



H = BN + H M + H Ze 



(1) 



where N is the molecular angular momentum excluding 
nuclear spins and if hi is the hyperfine Hamiltonian 24| . 



Hi 



hf 



i=l 



aN-Ii+ah-T-h+ali-h, (2) 



where Ii is the nuclear spin of nucleus i. The first term in 
the hyperfine Hamiltonian is the electric-quadrupole in- 
teraction, which represents the interaction of the nuclear 
quadrupole moment Q i with the electric field gradient V, 
at nucleus i. The second term is the nuclear spin- rotation 
interaction between the magnetic moment of the nuclear 
spin and the magnetic moment created by the molecular 
rotation. The third and fourth terms describe the tensor 
and scalar parts of the nuclear spin-spin interaction, re- 
spectively [25J. The interaction with the magnetic field 
is described by 



-g r fj, N N- B 
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gifx N Ii- B(l- ai), (3) 



which consists of rotational and nuclear contributions 
respectively, where g r and g^ are the rotational and 
nuclear (/-factors and <Xj is a nuclear screening factor. 
The constants c, have been calculated by Aldegunde et 
al. [24J. We use the fully uncoupled space- fixed basis, 
\I Rh M Rh )\Ic s M Cs )\NM N ), where /s7 Rb = f, / 133Cs = \ 
and the projections of the corresponding angular mo- 
menta are defined with respect to the z-axis directed 
along the magnetic field vector. 

To evaluate the outcome of collisions, we follow the 
approach described in Refs. [2?" 



27[. The total scatter- 



ing wavefunction, for molecules A and B, is expanded 
in the uncoupled space-fixed basis \taTbLMl), where 

\n) = \i£Im^JcI M i1 N ^ )M n)> for molecule i, and L 
is the end-over-end angular momentum of the molecule 
- molecule collision complex. This basis is then sym- 
metrized with respect to the interchange of two identical 



FIG. 1: (color online) The hyperfine structure of the 
87 Rb 133 Cs(X 1 E+) rotational ground-state. The hyperfine 
states considered here are highlighted in red. 



molecules |2j|[29j]. In a magnetic field, the total angular 
momentum projection M to t and the parity are conserved 
and the Hamiltonian can be block diagonalizcd. For each 
given M to t and parity, sets of coupled-channel equations 
are constructed, solutions to which are propagated out- 
wards over a radial grid of molecule - molecule separation, 
using the improved log-derivative method of Manolopou- 
los [30(. The radial grid was divided into short (4 to 
80 ao) and long range (80 to 1000 ao) regions, with grid 
step sizes of 0.0305 and 1 ao respectively. At 1000 ao, 
the log-derivative matrix is transformed into the basis 
of asymptotic eigenvectors and matched to the appro- 
priate scattering boundary conditions [3~fl ] to obtain the 
S'-matrix. The scattering cross sections are calculated 
from the S'-matrix elements (26j |. 

The dimension of the two-molecule basis set grows very 
fast as the number of the rotational states N and the or- 
bital angular momentum states L is increased. The ma- 
jority of the calculations presented here are performed 
with the basis set including three rotational states N < 2 
for each molecule and all L < 2 states. While this ba- 
sis set does not yield fully converged results, it is large 
enough to encapsulate the essential physics of molecule- 
molecule interactions [32| . In addition, we examined 
carefully the basis set dependence of our results. The 
calculations were carried out for M to t = ±10 and ±8, 
yielding 976 and 3129 coupled differential equations re- 
spectively. 

There is currently no full potential energy surface 
available for the RbCs(X 1 £+)-RbCs(X 1 £+) molecular 
system. We use the potential surface computed for 
NH( 3 £-)-NH( 3 Xr) collisions by Janssen et al. j^. In 
order to examine the potential surface dependence and 
present results relevant for the RbCs - RbCs system, we 
scale the NH - NH potential surface and compute the 
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FIG. 2: (color online) s-wave elastic and hyperfine relax- 
ation cross sections as a function of the collision energy. For 
two molecules in the initial (MRb = — §, Mas = — |) state the 
elastic (black-solid) and hyperfine relaxation (red-diamonds) 
cross sections are shown at a magnetic field of 10 Gauss. For 
two molecules in the initial (MRb = — ^,Mc s — — i) state, 
the relaxation cross section is shown for magnetic fields of 
10 (blue-circles) and 1000 (green-triangles) Gauss; the elastic 
cross section is very similar to the (M^h = — §,Mc s = — |) 
elastic cross section and is not shown. 



results in a wide rage of the scaling parameters. 

We being by analyzing the probability of nuclear spin 
relaxation of two molecules initially prepared in the maxi- 
mally spin-stretched state, labeled in Figure 1 as (Mjtb = 
— 3/2, Mcs = —7/2), and two molecules prepared in a 
lower energy state (Afiu, = — l/2,A/cs = —7/2). Figure 
2 shows the results of the scattering calculations, illus- 
trating two striking observations. First, the cross sections 
for nuclear spin relaxation have a large magnitude and 
can, for the initial state (A'/Rb = — 1/2, Mc s = —7/2) 
be comparable with the cross sections for elastic scat- 
tering. Second, the cross sections for the inelastic re- 
laxation from the (A/fib = —1/2, A/cs = —7/2) state un- 
dergo a rapid increase at a certain collision energy, mark- 
ing the opening of a barrierless collision channel. Below 
this collision energy all relaxation channels contain a d- 
wave barrier. This illustrates that there are significant 
couplings between the nuclear spin states. For collisions 
between molecules in the (A/Rb = — 3/2, Mq s = —7/2) 
state, the predominant coupling responsible for Feshbach 
resonances is the second order interplay of the nuclear 
quadrupole and potential interactions. For molecules in 
the (Mftb = — 1/2, Mcs = —7/2) state and in other non- 
strectched states an additional first order scalar spin-spin 
coupling also contributes. 

The same couplings give rise to magnetic Feshbach res- 
onances. In order to illustrate this, we compute the elas- 
tic and inelastic scattering cross sections for molecules 
prepared in the state (Afpib = — 1/2, Mc s = —7/2) as 
functions of the magnetic field. Figure 3 demonstrates 
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FIG. 3: (color online) Magnetic Feshbach resonances in 
molecule-molecule collisions computed using the N < 2 L < 2 
basis set at collision energy of InK. Upper panel: mag- 
netic Feshbach resonances for two molecules in low-field 
seeking (A/iu> = — s, AIc s = — \ ) state. Lower panel: two 
magnetic Feshbach resonances for two molecules in the 
(A^Rb = §,Mcs = |) absolute ground state, at 1210.1 and 
1670.2 Gauss. 



that the cross sections exhibit the resonant variation over 
the magnetic field interval as wide as ~ 10 Gauss. The in- 
elastic cross sections remain generally much smaller than 
the elastic cross section, as long as the collision energy 
is low enough (to the left of the enhancement point in 
Figure 2). 

The absolute ground state of RbCs, labeled in Figure 1 
as (MRb = 3/2, Mcs = 7/2) is particularly interesting for 
applications that need to avoid inelastic collisions. Fig- 
ure 3 illustrates that molecules prepared in this state ex- 
hibit Feshbach resonances that affect the scattering cross 
sections over an interval of 0.1 mGauss. For collisions 
between molecules in the absolute ground state, the cou- 
plings are the same as for the {Mub = — 3/2,A/c s = 
-7/2) state. 

The quantitative accuracy of the results presented is 
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FIG. 4: (color online) For two molecules in the absolute 
ground state (Mrtb = §,Mc s = I), the width of the largest 
Feshbach resonance in the range of magnetic fields from to 
4000 Gauss is shown as a function of the interaction potential. 



limited by the interaction potential surface and the basis 
set restrictions. While it is unlikely that increasing the 
basis set may decrease the resonance widths observed, 
it is necessary to confirm this. To perform calculations 
using larger basis sets, approximations are required to 
make the calculations feasible. The highest order off- 
diagonal coupling, that is predominantly responsible for 
Feshbach resonances, is the second order interplay of the 
nuclear quadrupole and potential interactions. To a first 
approximation we can neglect all states that are not cou- 
pled to a desired initial channel via this second order 
coupling. Practically this requires neglecting channels 
that differ from the initial channel by | AM/ tot | > 2iV max , 
where M hot = M$ + M%> + M<£> + and N max 

is the largest value of the rotational angular momentum 
included in the calculation. Neglecting these states ap- 
proximately halves the number of channels in the calcu- 
lation while only introducing a 1% error in the calculated 
cross sections. 

Using this approximation, we can incorporate addi- 
tional basis states. For molecules in the absolute ground 
state, increasing the number of partial waves from L < 2 
to L < 4 triples the number of Feshbach resonances to 
approximately three resonances per 1000 Gauss. While 
most of these additional resonances arc very narrow, 
there is always at least one resonance per 1000 Gauss 
with the width > 0.1 milliGauss. Increasing the number 
of rotational states from three to four (N < 3) rotational 
states per molecule, does not affect the resonance den- 
sity but increases the resonance widths by up to a factor 
of approximately thirty. Our predictions regarding the 
widths of the resonances should thus be considered as 
lower bounds. 

To demonstrate the effect of interaction potential vari- 
ation on the width of the resonances, we scale the NH- 



NH potential by a linear factor, V" scalcd = A sca i c V~ original . 
Figure [4] displays the width of the widest resonance for 
molecules in the absolute ground state in the magnetic 
field range between and 4000 Gauss as a function of 
the potential scaling factor, using the above approxima- 
tion with three rotational states per molecule and L < 2. 
The resonance width oscillates between 10~ 7 Gauss and 
10~ 2 Gauss, with an average of 3.6 x 10~ 4 Gauss. As 
the potential is varied the larger resonances are moved 
out of the magnetic field region considered, and thus for 
the smaller values in figure H] a larger resonance could be 
found by considering a larger range of magnetic field. 

In summary, we have presented a rigorous analysis of 
the effects of magnetic fields on collisions of molecules. 
We showed that the interaction of the nuclear spins with 
the magnetic field gives rise to magnetic Feshbach res- 
onances mediated by an interplay of hypcrfine interac- 
tions and the electronic interaction potential between 
molecules. Our calculations show that molecules pre- 
pared in the absolute ground state may exhibit mag- 
netic Feshbach resonances with widths > 1 milliGauss. 
Molecules prepared in excited Zeeman levels of the ro- 
vibrationally ground state may exhibit even broader res- 
onances, with widths > 1 Gauss. Although excited Zee- 
man levels are subject to inelastic Zeeman relaxation 
changing the nuclear spin projections, we found that cer- 
tain states exhibit broad Feshbach resonances and sup- 
pressed inelastic relaxation. These states appear to be 
ideal for experiments exploiting the magnetic field control 
of intermolecular interactions. These states exhibit the 
unusual dynamics of inelastic collisions that is extremely 
sensitive to the opening of excited collision thresholds, 
which can be used as a new mechanism for controlling 
binary molecule interactions. 

This work was supported by NSERC of Canada 
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